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Summary 

Reaction-rate coefficients and thermodynamic and transport properties are reviewed and sup- 
plemented for an 11-species air model. These coefficients and properties can be used for analyzing 
flows in chemical and thermal nonequilibrium up to temperatures of 30000 K. Such flows will likely 
occur around currently planned and future hypersonic vehicles. Guidelines for determining the state 
of the surrounding environment are provided. Curve fits are given for the various species properties 
for their efficient computation in flow-field codes. Approximate and more exact formulas are pro- 
vided for computing the properties of partially ionized air mixtures in a high-energy environment. 
Limitations of the approximate mixing laws arc discussed for a mixture of ionized species. An elec- 
tron number-density correction for the transport properties of the charged species is given. This 
correction has generally been ignored in the aerospace literature. 

Tabulated values of the curve-fit coefficients and computer subroutines to evaluate the various 
species properties using these coefficients arc available from the NASA Computer Software Manage- 
ment and Information Center (COSMIC). 

Introduction 

Currently envisaged transatmospheric and aeroassist missions (refs. 1 to 4) have created a 
resurgence of interest in the aerothermodynamic design of hypersonic vehicles. However, the 
velocities and altitudes at which these proposed craft would operate are different, and sometimes 
more severe, than those experienced in the past. As a result, the nonequilibrium flow environment 
that will surround these vehicles will considerably impact the vehicle aerodynamics, thermal loads, 
and propulsion-system efficiency. Since such an environment is difficult to simulate in current ground- 
based test facilities, the design of these future vehicles will rely heavily on numeric calculations. In 
turn, these calculations will require a good understanding of the physical modelling required to 
simulate these phenomena. 

Under hypersonic flight conditions, a vehicle travelling through the atmosphere will excite the air 
that flows around the body to very high temperatures as the kinetic energy of the vehicle is dissipated 
to the gas. Depending on the flight velocity, various chemical reactions will be produced behind a 
shock wave as shown in figure 1 (which is adapted from ref. 5) for the stagnation region of a sphere 
with a radius of 30.5 cm. These reactions will affect the properties of air and cause considerable 
deviation from those of a thermally and calorifically perfect gas. A vehicle flying through the higher 
reaches of the atmosphere at high velocities may also experience thermal nonequilibrium (fig. 1), 
since the lower density reduces the collision frequency, and the high velocity results in smaller 
transit times for the air molecules. Both of these processes create a delay in the equilibration 
of translational, rotational, vibrational, and electronic modes of the thermal energy. Under these 
conditions, the modelling of the air chemistry requires a multitemperature approach in contrast to 
classical single-temperature formulations. 

Four regions (I to IV) are delineated in figure 1. These regions show when the various chemical 
activities are initiated at a given altitude and velocity. The initiation of chemical and thermal 
nonequilibrium processes for different velocity and altitude conditions is similarly depicted through 
regions A, B, and C. This figure clearly shows that the set of chemical reactions and thermodynamic 
and transport properties would change continuously for a given flight trajectory. For example, in 
regions A and B (i.e., before initiation of thermal nonequilibrium), the specific heat at constant 
pressure C p would change as follows: 

C p = Constant in region I before the excitation of vibrational energy mode 

C p = C P (T) in region I after the excitation of vibrational energy mode and 
before the dissociation of oxygen 

C p = C p {Ci,T) after the dissociation of oxygen 


Similarly, the equation of state changes along the flight trajectory as the thermal equilibrium and 
thermal nonequilibrium regions (ref. 6) are traversed: 

p = T (in regions A and B) 

M 

^univ rrt (' * r*\ 

P = T tr ( m re S lon C) 

where 7t r is the translational temperature in an environment of thermal nonequilibrium. 

In numerical simulations, the thermodynamic and transport properties and reaction-rate coef- 
ficients (in the case of finite-rate chemistry) are usually required. It is obvious from the previous 
discussion that these properties and the equation of state should be evaluated carefully when chemi- 
cal and thermal nonequilibrium conditions exist in the flow field around a hypersonic vehicle. Under 
chemical and thermal equilibrium conditions, the transport and thermodynamic properties of high- 
temperature air and its components are well documented in the literature (refs. 7 to 12). However, 
for flows with finite-rate chemistry, the individual species properties and appropriate mixing laws 
that are required are not as well established. For example, in a partially ionized gas mixture, the 
conventional mixing laws (refs. 13 and 14) developed for nonionized mixtures cannot be extended to 
higher temperatures without considerable error (ref. 15). 

The purpose of this report is to review the thermodynamic and transport properties and the 
reaction-rate coefficients of the most important reactions for the 11 constituent species of air (N, O, 
N 2 , O 2 , NO, N + , 0 + , N^, Oj , NO + , e“) for temperatures up to 30000 K. Those species properties 
that are not available in the literature for this 11-spccies model are provided, and curve fits are 
obtained for all properties to permit their efficient computation in flow-field codes. Tabulated values 
of the curve-fit coefficients and computer subroutines to evaluate the various species properties using 
these coefficients are provided in reference 16. These values and codes may be obtained for a fee 
from: 

COSMIC 

Computer Services Annex 
University of Georgia 
Athens, GA 30602 
(404) 542-3265 


Request the code by the designation LAR 14447. This code is written in FORTRAN 77 for use on 
computer with a FORTRAN compiler. 

Approximate and more exact mixing laws for the various species properties are also provided 
for partially ionized gas mixtures. The limitations of the approximate mixing laws are pointed out, 
especially for a mixture of ionized species. An electron number-density correction for the transport 
properties of the ionized species, which has been generally neglected in the aerospace literature, 
is provided. Sources of the input data used in the calculation of various flow-field properties are 
identified. Appropriate formulas are provided for using these properties in computations of flows 
with thermal nonequilibrium. 
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curve-fit coefficients for collision cross section ratio B*j 
(eq. (47)) 

coefficient in modified Arrhenius form of backward 
reaction-rate constant (eq. (3b)) 
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curve-fit coefficients for diffusion coefficient D t j 
(eq. (42c)) 

coefficient in modified Arrhenius form of forward 
reaction-rate constant (eq. (3a)) 

matrix elements of first Chapman-Enskog formula 
ratio of collision cross sections, 

j l J 

curve-fit coefficients for equilibrium constants K cqT 
(cqs. (5c) and (5d)) 

curve-fit coefficients for frozen thermal conductivity of 
species i,Kji (eq. (26)) 

coefficients of polynomial curve fits for thermodynamic 
properties, n = 1,2, 7 (cqs. (10a) to (10c)) 

curve- fit coefficients for viscosity of species i, //, 

(cq. (25)) 

curvc-fit coefficients for collision cross sections ft' ' ^ 
(eq. (45)) 
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curve-fit coefficients for collision cross sections £l\-' 

(eq. (46)) 

first Bohr radius, 0.52918 x 10 -8 cm 

temperature exponent for backward reaction-rate 
constant 

temperature exponent for forward react ion- rate constant 


ratio of collision cross sections, (5fl|j^ — 4n|j’ 3 ^) j 


(bl) 


mass diffusion parameter defined by equation (A 10) 
specific heat at constant pressure, (§t) p ’ cal/g-molo-K 
reactive specific heat at constant pressure, 

, cal/g-molc-K 


frozen specific heat at constant pressure, 

, cal/g-mole-K 

p 

specific heat at constant pressure of species i, 

(fr) 1 cal/g-mole-K 

electronic component of C p i , cal/g-mole-K 
internal component of C p ^ , cal/g-mole-K 
rotational component of C pd , cal/g-mole-K 
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Cl 

Fi 
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K e 

Fe\ 


K, 


eq,r 


K f 


translational component of C PsZ -, cal/g-mole-K 
vibrational component of C p<? ;, cal/g-mole-K 

specific heat at constant volume, i cal/g-mole-K 

specific heat at constant volume of species z, 

(fr) ’ cal/g-mole-K 

internal component of C v ^ cal/g-mole-K 
translational component of C v i, cal/g-mole-K 
coefficient of self diffusion, cm 2 /sec 
binary diffusion coefficient, cm 2 /sec 
= pD{j , cm 2 -atm/scc 

multicomponent diffusion coefficient, cm 2 /sec 

activation energy for backward reaction r, erg/g-mole 

activation energy for forward reaction r, erg/g-mole 

internal energy of mixture, cal/g-mole 

internal energy of species z, cal/g-mole 

free energy of species z, cal/g-mole 

free energy of species i at the standard state of 1 atm 
pressure (1 atm — 101.3 kPa), cal/g-molc 

NS 

enthalpy of mixture ^2 Cih z , cal/g-mole 
i= 1 

Planck constant, 6.6261 x 10~ 27 erg-sec 

enthalpy of species z, cal/g-mole 

specific-heat function for species z, cal/g-mole-K 

standard heat of formation of species i at temperature 
T re f, cal/g-mole 

conversion factor from calories to ergs, 4.184 x 10 7 
ergs / calorie 

fcth component of diffusion mass flux of species z, k = 1 
to 3, g/cm 2 -sec 

total effective thermal conductivity of mixture in 
thermodynamic equilibrium, Kj + K r , cal/cm-sec-K 

thermal conductivity of electrons, cal/cm-sec-K 

component of thermal conductivity of mixture due to 
electron excitation, cal/cm-sec-K 

equilibrium constant for reaction r, kj^/k^^ 

frozen thermal conductivity of mixture in thermo- 
dynamic equilibrium, K ir + K mi , cal/cm-sec-K 
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K f,i 

K'mt 
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K r 
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frozen thermal conductivity of species i in thermo- 
dynamic equilibrium, K trp + cal/cm-sec-K 

internal component of frozen thermal conductivity of a 
mixture in thermodynamic equilibrium, K rot + K v fo + I<eh 
cal/cm-sec-K 

internal component of frozen thermal conductivity of 
species i in thermodynamic equilibrium, cal/cm-sec-K 

reaction component of thermal conductivity of a mixture 
(eqs. (A19) and (A20)), cal/cm-sec-K 

translational thermal conductivity of mixture from first 
Chapman-Enskog approximation, cal/cm-sec-K 

translational thermal conductivity of mixture without 
contributions caused by electron-heavy particle collisions, 
K\ t — K c , cal/cm-sec-K 

translational component of thermal conductivity of 
species i> cal/cm-sec-K 

vibrational component of thermal conductivity of mix- 
ture, cal/cm-sec-K 

Boltzmann’s constant, 1.38066 x 10“ 16 crg/K 

backward reaction-rate coefficient for reaction r, 
cm 3 /mole-sec or cm 6 /mole 2 -sec 

forward reaction-rate coefficient for reaction r, 
cm 3 /mole-sec 

Lewis number, defined by equation (A9) 
frozen binary Lewis number, pC P jDij/ K j 
total binary Lewis number, pC p D t j/K 
multicomponent Lewis number, pCpjDjj j K j 
molecular weight of mixture, g/g-mole 
molecular weight of species z, g/g-mole 
mass of particle z, g 

Avogadro’s number, 6.0221 x 10 23 molecules/g-mole 

total Prandtl number, C p p/K 

frozen Prandtl number, C P jp/Kj 

number density, particles/cm 3 

number density of species z, particles/cm 3 


P 


pressure, atm 


Pe 


electron pressure, atm 


Pe m 


limiting value of electron pressure, atm 
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translational partition function per unit volume, 1/cm 3 

fctli component of overall heat-flux vector, k — 1 to 3, 
cal/cm 2 -sec 

universal gas constant, 1.987 cal/g-mole-K 

radius of electron orbit from solution of Schrodinger 
equation, cm 

Coulomb cutoff (or Debye shielding) radius, cm 
density cutoff radius, cm 

entropy of species i at temperature T, cal/g-mole-K 

temperature under thermodynamic equilibrium, K 

characteristic reaction temperature for backward reac- 
tion r, Ebf./k, K 

characteristic reaction temperature for forward reaction 
r, Ej r /k, K 

electron temperature, K 

electronic excitation temperature, K 

reference temperature, 298.15 K 

rotational temperature, K 

translational temperature, K 

vibrational-elcctron-electronic excitation temperature, K 
vibrational temperature, K 
time, sec 

mean molecular velocity of species z, cm/sec 

A*th component of diffusion velocity of species z, k = 1 
to 3, cm/sec 

mass rate of formation of species i, g/cm 3 -sec 
concentration of species z, moles /volume 
mole fraction of species i 

A th component of general orthogonal coordinate system, 
k = 1 to 3, cm 

constants (used in eq. (2e)) given in table I 
stoichiometric coefficient for reactant i in reaction r 
stoichiometric coefficient for product i in reaction r 
ratio of specific heats, C p /C v 
mole mass ratio of species j, Xj/p 

defined by equation (34), cm-sec 
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Subscripts: 

6 

e 

el 

i 

3 

l 

P 

r 

rot 

tr 

ve 

1/ 

Abbreviations: 

ASTV 

mol. 

NASP 

NATA 


defined by equation (35), cm-sec 

mean- free path of pure species i, cm 

modified mean-free path of species i in a mixture, cm 

mean-free path of species i in a mixture, cm 

viscosity of mixture from first Chapman-Enskog approxi- 
mation, g/cm-sec 

viscosity of species z, g/cm-sec 

density of mixture, g/cm 3 

density of species z, g/cm 3 

average collision cross section (used for diffusion, viscos- 
ity, translational, internal, and reaction components of 
thermal conductivity) for collisions between species i and 
j. A 2 (1 A= 10" 8 cm) 

average collision cross section (used for viscosity and 
translational component of thermal conductivity) for 
collisions between species z and j, A 2 

average collision cross sections (used for translational 
component of thermal conductivity) for collisions be- 
tween species i and jf, A 2 

backward reaction 
electron 

electronic excitation 
species i 
species j 
species l 

constant pressure 
reaction 

rotational energy mode 
translational energy mode 
vibrational-electron-electronic energy mode 
= NS, total number of species 

Aeroassisted Space Transfer Vehicle 
molecules 

National Aero-Space Plane 
Nonequilibrium Arc Tunnel Analysis 



NIR 


total number of independent reactions 

NJ sum of reacting species (NS) plus number of catalytic 

bodies 

NR total number of reactions 

NS total number of species 

Chemical Kinetic Model and Reaction Rates 


Thermal Equilibrium 

When chemical reactions proceed at a finite rate, the rate-of-production terms appear h{ the 
energy equation, when formulated in terms of temperature, and in the species continuity equations 
(refs. 17 to 19). For a multicomponent gas with NS reacting chemical species and NR chemical 
reactions, the stoichiometric relations for the overall change from reactants to products are 


NJ k f' r NJ 


where r = 1, 2, NR and NJ is equal to the sum of the reacting species (NS) plus the number 
of catalytic bodies. The quantities aj >r and /?^ r are the stoichiometric coefficients for reactants and 
products, respectively, and kj r and k^ r are the forward and backward rate constants. The quantities 
X x denote the concentrations of the chemical species and catalytic bodies in moles per unit volume. 
The catalytic bodies (NJ - NS) may be chemical species or linear combinations of species that do 
not undergo a chemical change during the reaction. 

The net mass rate of production of the zth species per unit volume resulting from all the reactions 
NR may be obtained (ref. 18) from 


Wj = M { 


NR / j y \ NR / j x 

5 (£),-£ (£), 


or 


NR 

Wi — Afj y ^ (fli, r ~~ a i y r)( R /,r — R b,r) 
r=l 


where 


NJ 


R f r = */,r II (Wf* 

j =1 


NJ 


R b,r = hr II 

J = l 

Here, the mole-mass ratio 7^ (also known as the mole number) is defined as 


£ = j| O' = 1,2,..., NS) 

7 3 = \ NS 

H U-w + 1 NJ) 

0"NS).« 


(2a) 


(2b) 


(2c) 

(2d) 


(2e) 
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The constants are functions of the catalytic efficiencies of the NS species and are 

determined from the linear dependence of the catalytic bodies upon the NS species. Values of these 
constants for the 11-species air model are given in table L 

The reaction rates in equation (1) or equations (2c) and (2d) are expressed in the modified 
Arrhenius form as 


k f,r = A f.rT B f’ r exp(-T Dfr /T), 
h,r = \rT BbT exp (~T Dbr /T), 



(3a) 

(3b) 


where 


NJ 

a r = - 1 

i = 1 

NJ 

Pr = Pis ~ ^ 

i = 1 


(4a) 

(4b) 


and where Tdj t an< 3 ^D br are the characteristic reaction temperatures for the forward and backward 
reactions, respectively. Values for the reaction rates and k^ r are tabulated in table II for the 
11-species air model. For a specified temperature, density, and species composition, equations (2) to 
(4) can be used to obtain the production rate of a species i in a multicomponent gas by employing 
the catalytic body efficiencies and reaction rates from tables I and II. The first seven reactions and 
reaction rates in table II are taken from reference 20 and were employed in reference 17 for the 
7-species air model (N, 0, N2, O2, NO, NO + , e“). Reaction rates for reactions 8 to 20 of table II are 
taken from reference 21. Some of these reactions have been regrouped here (and in ref. 22) through 
the use of third bodies M\ to A/4, which is similar to the approach used in reference 17. 

The reaction rates given in table II are appropriate for flow velocities of about 8 km/sec (i.e., up 
to the Shuttle-type reentry velocities). For higher flow velocities, the backward reaction rates can 
be obtained from the forward reaction rates 1 by using the following relation: 


— 


kf t r(T) 

K eq ,r(T) 


(5a) 


The equilibrium constants K eq , r have been obtained by using the atomic partition functions and 
the molecular partition functions provided in reference 23. The electronic partition function for the 
atomic species and internal partition functions of molecular species are tabulated in reference 23. The 
internal partition functions of the molecular species were obtained using existing data on molecular 
constants augmented by solving the Sc.hrodinger equation. The translational component required 
for obtaining the total partition function (which is used to compute the equilibrium constant) is 
obtained from 

Qti,i = (2tt m t kT/h 2 f/ 2 (5b) 

where Qtr,z ls the translational partition function per unit volume and h is the Planck constant. 

The computed values of K e q,r have been curve- fit ted here by the least-squares curve-fit method 
as a function of temperature using the expression 


K eq , r — | exp (Ftf cqr ) j (10 A /T)^ Ak ^ tZ + ^eq.r^ 3+ ^e qi r^ 2+ ^^q,r Z+£: /vcq,r) (5 C ) 


or, 


ln(*W) = A KoqT Z 5 + + C Kcfx Z* + D h - cqr Z 2 + E Kcq Z + F Kcqr 


(5d) 


1 The high-temperature data for the forward reaction rates are not well-defined. 
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where 


Z — In (10 4 /T) (5c) 

The curve-fit coefficients in equations (5c) and (5d) are given in table III. Caution should be exercised 
in evaluating the equilibrium constant K eqr from the coefficients given in table III. All the six 
significant digits for the coefficients should be used. Equation (5d) is preferred to equation (5c) in 
computing the equilibrium constant for reasons of accuracy, because the exponent of the exponent ial 
term in equation (5c) is a large number. Since the electronic partition functions of atomic species 
(tabulated in ref. 23) are significantly affected at high temperatures and low densities, the curve-fit 
coefficients are given for 6 different values of the total number density. These number densities 
cover the range of practical interest for aerospace applications. Figure 2(a) shows the variation 
of the equilibrium constant with temperature for different values of the total number density. The 
equilibrium constant clearly exhibits a number-density dependence at very high temperatures and low 
total number densities. Figure 2(b) gives a comparison of the present curve fit with the exact values 
calculated from the partition functions. Values obtained from Park's curve fits (ref. 23) and from 
table II by employing the relation K v q , r = kf, r /kbj- are also shown. Park’s five curve-fit coefficients 
were obtained by using five discrete temperatures (2000, 4000, 6000, 8000, and 10 000 K); thus, they 
are in disagreement with the exact values and the present curve fits at the higher temperatures. 
Comparison of the exact values with those obtained from table II also gives an indication of the 
temperature range beyond which the backward reaction rates of table II should not be employed. 

The equilibrium constant /v e q,r is usually given as a function of temperature only. Figure 2(a) 
and table III indicate, however, that under low-density and high-temperature conditions K eqr is also 
a function of the total number density n. This dependence arises because the electronic partition 
functions of nitrogen and oxygen atoms and ions are dependent upon the total number density under 
those conditions. In reality, there is a finite (as opposed to infinite) number of electronic states that 
contribute to the electronic partition Function. This finite number of states is obtained because 
(a) the orbits of electrons cannot extend beyond the distance to the nearest neighboring particle 
(a phenomenon known as the density cutoff) and (b) the Coulomb field is truncated due to the 
perturbation of the field by the electrons and ions in the neighborhood (a phenomenon known as 
the charge shielding). On the average, the distance to the nearest particle is given in centimeters as 

r d = 1/n 1 / 3 (5f) 

where n is the total number density in particles/cm'* and r ( [ is the radius of an electron orbit as 
the result of density cutoff. In the hypersonic flight regime of interest, this radius prevails over the 
Coulomb cutoff radius r c , (or Debye shielding) defined in centimeters by 

f T \ 1/2 

» - 6.90 ( — J (5g) 

Further, the quantum theory gives the allowable radius r a of an electron orbit as follows: 

r« = a 0 N 2 (5h) 


where aq = 0.52918 x 10 -8 cm and is known as the first Bohr radius, and N (which is an eigenvalue 
in the solution of the Schrodinger equation) is called the principal quantum number. Equating the 
smaller of the two radii given by equations (5f) and (5g) with equation (5h) yields the cutoff quantum 


number N m ax as 


A max 



(5i) 


Only electronic states with principal quantum numbers up through A r max contribute to the electronic 
partition function. Since the contribution from the highly excited states with quantum numbers near 
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Nmax is proportional to A^* ax exp(— E/kT), where E is the ionization energy, these highly excited 
states dominate the overall partition function at sufficiently high temperatures or low total number 
densities, and the partition function then becomes strongly dependent on number density (eqs. (5f), 
(5g), and (5i)). Under these conditions, most of the atoms present in the gas are in the highly excited 
electronic states near N mRX rather than in the ground state. Since these highly excited atoms have 
very different properties than ground-state atoms, the reaction-rate data given in table IT, which were 
determined for primarily ground-state atoms, are no longer appropriate for the low-density, high- 
temperature conditions at which the number- density dependence of the partition function becomes 
significant. Therefore, more general models, which account explicitly for the behavior of the excited 
atoms, are required. This regime lies outside the range of primary interest for aerospace applications, 
however, and is not covered in this study. 

Park (ref. 24) has recently presented a new set of reaction rates for use in comput ing highly ionized 
flow in the hypersonic environment. In particular, the power of the preexponential temperature of the 
rate coefficients for the associative ionization and dissociative recombination reactions (reactions 7, 
10, and 13 in table II) is assigned a value of zero at temperatures greater than 6000 K to keep 
the rate coefficients from becoming unrealistically large at very high temperatures, for which no 
experimental data are available. Since the power of the preexponential temperature of the rate 
coefficients is negative and zero for reactions 10 and 13, respectively, in table II, Park's suggestion 
is already implemented in the reaction rates for these reactions. The preexponential power of the 
rate coefficient for reaction 7 may need to be reexamined when the temperatures become excessively 
large. Further, for extremely high temperatures, such as those encountered with earth reentry 
velocities greater than 12 km/sec, much slower reaction rates have been recommended (refs. 25 
and 26) for electron-impact ionization reactions 8 and 9 given in table II. Also, the reaction rates 
given in reactions 8 and 9 in table II are from expansion flow data, which tend to be lower than data 
obtained under compressive flow conditions. 

Thermal Nonequilibrium 

The reaction rates given in table II were originally used by Blottner (ref. 17) and Dunn and 
Kang (ref. 21) in the context of a single temperature assuming thermal equilibrium. Park’s 
(ref. 25) guidelines (in the context of his two-temperature model) may be used for defining the rate- 
controlling temperature in dissociation and the electron-impact ionization reactions under thermal 
nonequilibrium conditions. These guidelines were used in reference 27, for example. 

Based on the preferential dissociation concept, Park has suggested the use of a temperature 
weighted with the vibrational temperature to characterize dissociation reactions (ref. 25). The 
reaction rates in Park’s model are assumed to be dictated by the geometric average temperature 

= x/Ttr^vib (6a) 

and the dissociation reaction rates are given by 

= ^/,r^av^ CXp(-T D/r /r av ) (6b) 

The recombination (or backward) reaction rate, however, depends only on the temperature of the 
impacting particles and may, therefore, be evaluated from 

kb ,r ~ kf r (T iT )/h eqT (T tr ) (6c) 

Treanor and Marrone (refs. 28 and 29) suggested a more rational (but slightly more difficult) way 
than Park s method to account for the effect of vibrational relaxation on dissociative reactions with 
the preferential dissociation concept. They suggested the use of a vibrational coupling factor (refs. 27 
to 29) with the dissociation reaction rates obtained under the assumption of thermal equilibrium. 
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Recent work of Jaffe (ref. 30), based on collision theory and using methods of statistical 
mechanics, yielded no evidence of preferential weighting to any particular energy mode in obtaining 
the total energy available in a collision, whether it is an elastic, inelastic, or reactive encounter. Jaffe 
found that the multitemperature effects on the reaction rates were small for dissociation. These 
findings were supported by those of Moss et al. (ref. 31), who carried out flow-field analyses with 
the Direct Simulation Monte Carlo (DSMC) approach. Thus, a weaker dependence of kj r on T vib 
(such as the one suggested by Sharma et al. (ref. 32) as T av = T^r^J) mi ght be more realistic, 
especially for highly energetic flows. 

It is obvious that the multitemperature kinetic models for high-energy flows based on both 
preferential and nonpreferential dissociation assumptions employ a considerable degree of empiricism. 
They exemplify the degree of uncertainty that exists in modelling the multitemperature kinetics. 
Quantum mechanical studies of the type in reference 30, supplemented by nonobtrusive laser 
diagnostic studies, would be desirable for establishing these models on a sounder basis. 

Species Thermodynamic Properties and Mixture Formulas 

Thermal Equilibrium 

Thermodynamic properties (i.e., C p j and h z ) are required for each species considered in a finite- 
rate flow-field calculation. For calculations with chemical equilibrium, the free energies Fj are 
also required. Since the multicomponent gas mixtures are considered to be mixtures of thermally 
perfect gases, the thermodynamic properties for each species are calculated by using the local static 
temperature. Then, properties for the gas mixture are determined in terms of the individual species 
properties through the following relations: 


NS 

h — ^ ^ Cfoi 
1 


(7a) 


where 


and 


where 




(7b) 


(8a) 

(8b) 


The mixture C Pj as defined by equation (8a) is the “frozen” specific heat. This definition does 
not account for species production or conversion due to chemical reactions. Frozen specific heat 
is commonly used in defining the Prandtl and Lewis numbers 2 for a mixture and is related to the 
mixture enthalpy h through the relation 


C Pf = 





i=l 



p 


NS 

Eocw 

i=l 


( 9 ) 


Expressions for C p j using the partition-function approach were obtained in reference 33, whereas 
a virial coefficient method was used in references 34 and 35. If the effects of nonrigidity of the rotor 


2 See appendix A for various definitions. 
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and anharmonicity of the oscillator are included in these two methods, the two formulations are 
equivalent (see section 3.1 of ref. 36) and can be used interchangeably. Using the virial formulation, 
Browne (refs. 34 and 35) obtained the thermodynamic properties as corrections to those of the 
monatomic gas in terms of the first and second virial coefficients and their temperature derivatives. 

Reference 33 provides curve fits for C p i and whereas references 34 and 35 provide tabulations 
of these data. Since the use of curve fits reduces the expense of computing the original functional 
relations, the tabulated thermodynamic properties of references 34 and 35 have been curve-fitted 
here as a function of temperature for the temperature range of 300 K < T < 30000 K. The following 
polynomial equations are employed for these curve fits: 

Specific heat: 

= A { + A 2 T + A 3 T 2 + A 4 T 3 + A b T 4 (10a) 

*Miniv 


Specific enthalpy: 


hi 4 , a 2 t , A 3 T 2 , a 4 t 3 , a 5 t 4 
R—f = A ' + ^r + ~ + — + — 



For equilibrium calculations, the following curve fit for the free energies F t may be used: 


(10b) 


^nm\T 


= A 1 [l-\n(T)} 


a 2 t a 3 t 2 a 4 t 3 

~2 6 12 ~ 


a 5 t 4 

20 



(10c) 


where is the free energy of species i at a pressure of 1 atin (standard state). The specific-heat 
data are easily curve-fitted to the polynomial form of equation (10a). Other polynomial forms are 
based on the following thermodynamic relations at constant pressure: 


dh{ — C Pi dT 

(lla) 

dSi = c Pl p 

(lib) 

dhi — T dS t 

(11c) 

dFi = — Si dT 

(lid) 


where F x is the free energy and Sj is the entropy of species i at temperature T. 

The polynomial coefficients have been evaluated by using the least-squares curve- fit technique. In 
particular, the following polynomial for the specific enthalpy has been curve-fitted to the tabulated 
values of references 34 and 35: 

h,--^ 6 fl un j v =] 4 1+ x 2 T + A 3 T 2 + A 4 T 3 + A 5 T 4 (12a) 

Constants A\ to A§ in equations (10a) to (10c) are then obtained from 

A n = nA n {n = 1,2, 3, 4, 5) (12b) 

The constant Aq in equations (10b) and (10c) is computed separately. This constant is related 
to the heat of formation through the relation 

^6^univ = i h i)r= 0 = r rcf - [(hj) Tref - (h,) T=0 ] (12c) 

In equation (12c), [(^i)r ref ~ (^j)r=o] available from the tabulated data to be fitted. Sources giving 
tabulated values of specific enthalpy and free energy generally use 0 K as the reference temperature. 
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To transform this reference temperature to the reference of 298.15 K employed here, it is necessary 
to know the heat of reaction at the new reference. Appendix B gives the method used for calculating 

(A/(/)r =2 98.15- 

It is sometimes useful to define a specific-heat function /q such that 


fk = 


hi - (Ahj )r r( , f 


( 1 2d) 


This relation may be used to obtain T from h, in an iterative manner. Equation ( 1 2d) may be 
expressed in a polynomial form by using equations (12a) and (12c) as follows: 


hi = (Aj + A 2 T + A 3 T :i + .4jT 3 + A b T q )R univ - 




(12e) 


T 

or, with the approximation f 0 rvf C p j dT ~ C p jT xv f. 


A\r> Cp/I 


hi = (A\ + A 2 T + A^T + A,\T 6 + A 5 T^)R nuW - 


T 


(12f) 


The constant Aj has been obtained by subtracting the remaining terms of the free-energy 
polynomial (eq. (10c)), evaluated with the other known constants, from the tabulated frce-cnergy 
data as follows: 


A 7 = A i[l — ln(T)] -^-T 


A 2rr A 3 T 2 A 4 T 3 


12 


A 5 T ’ 4 
20 


f? ~ (hj ) T = c 

«univT 


(13) 


Table IV is a listing of the polynomial constants {A\ to A-j) for the 11-specics air model. The 
polynomial coefficients have been obtained for five temperature ranges between 300 and 30000 K. 
To assure a smooth variation of thermodynamic properties over the entire temperature range, values 
of A\ to A 7 should be linearly averaged across the curve-fit boundaries (i.e., 800 < T < 1200, 5500 < 
T < 6500,14 500 < T < 15 500, and 24 500 < T < 25 500), because the curves are not continuous 
at the boundaries. A sample subroutine that evaluates the polynomial curve fits and performs the 
linear averaging is presented in appendix C. This routine may be easily modified to suit the user’s 
requirements. 

Temperature T in equations (10) to (12) is in kelvins. With the universal gas constant in cal/g- 
molc-K, the specific heat and enthalpies have the units of cal/g-mole-K and cal/g-mole, respectively. 

Figure 3 is a comparison of the values of specific heat obtained from the polynomial curve fit 
(eqs. (10)) with the data of Browne (ref. 35). Values provided by Hansen (ref. 37) are also included 
for comparison. Hansen’s values begin to deviate from those of references 33 and 35 beyond 4000 K 
for O2 (fig- 3(a)), beyond 8000 K for N2 (fig. 3(b)), and beyond 6000 K for NO (fig. 3(c)). This 
deviation may be a result of the rigid-rotor and harmonic-oscillator partition function employed by 
Hansen and of the neglect of electronic cont ribut ion. At higher temperatures, the population of levels 
corresponding to the nonparabolic regions of the potential-energy curve arc no longer negligible, and 
it becomes necessary to introduce the nonrigidity and anharmonicity corrections into the energy 
levels of the molecules as was done in reference 35. ,. 

Reference 38 has provided tabulated values for specific heats, enthalpies, and internal energies 
for the 11-spccics air model up to temperatures of 6000 K. These values arc nearly identical to those 
given in references 33 to 35. The minor differences that do exist may be due to different forms of the 
partition function, different spectroscopic data, inclusion of excited-state data, inclusion of isotopic 
effects, different heats of formation, or other factors. 
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Thermal Nonequilibrium 

The thermodynamic property data of references 33 to 35 can be used to separate the cont ributions 
of different internal energy modes to the specific heat, as is required in the case of multitemperature 
flow- field models. For example, in a two- temperature model (refs. 25 and 27), advantage may be 
taken of the fact that the translational and rotational energy modes are fully excited and equilibrated 
at room temperature; therefore, the heat capacities for these modes are independent of temperature. 
The combined vibrational-electronic specific heat for species (C pi ) yQ , can then be evaluated by 
using the value for the total specific heat C pp evaluated at temperature T ve and by subtracting out 
the constant contribution from the translational and rotational specific heats. This can be described 
(see fig. 4 for the various contributions) by 


[(CW)ve h-c = (Cp.iJlW, - (Cp,)tr - (C Pii ) rot (14) 

where the translational component of the specific heat (C p / )i r is | i? mi j v for all species, while the 
rotational component (Cy) r ot I s i?univ for diatomic species and 0 for monatomic species. 

The enthalpy h r for a two- temperature model can be evaluated similarly, since contributions 
from the translational and rotational modes are linear with temperature. Therefore, the vibrational- 
electronic enthalpy for species i, (hj) v c , can be obtained from the specific enthalpy hj evaluated at 
temperature T vc by subtracting the contribution from the translational and rotational enthalpies 
evaluated at T ve as well as the enthalpy of formation as follows: 


[(^)vc]r vr “ [^i] r vc [(C/4,/)tr + (Cp,i) rot](7\ r e — ^ref) — 


(15) 


The specific enthalpy from all the contributions of internal energy modes can then be obtained by 
adding the contributions of translational and rotational enthalpies (evaluated at the translational- 
rotational temperature) and the enthalpy of formation to the vibrational-electronic enthalpy as 
follows: 

hi(T, T ve ) = [(hi)y c]r w + {(Cp,ih + (C p ,i)rot}(T - T rcf ) + (A h{ )r n , f (16) 

The procedure outlined here for obtaining the specific heats and enthalpies for different internal 
energy modes can be used only with one- or two- temperature models. This approach may be used 
for a three-temperature model only if one of the energy modes is partially excited and the rest are 
fully excited. If more than one internal energy mode is only partially excited, the partition-function 
approach (with appropriate corrections for the rotor nonrigidity and oscillator anharmonicity) would 
be needed to obtain the thermodynamic properties of different energy modes. 


Species Transport Properties and Mixture Formulas 


Thermal Equilibrium 

The transport properties required in flow-field calculations are viscosity, thermal conductivity, 
and diffusion coefficients. The collision cross sections required for these properties have been 
recomputed herein using the same molecular data used previously by Yos (refs. 7 to 9). The 
computational techniques employed in the calculations are described in references 39 to 41; these 
references give details of the NATA (Nonequilibrium Arc Tunnel Analysis) code. In NATA, the 

average collision cross sections n for the collisions between species i and j are calculated from 
basic cross-section data as functions of temperature and gas composition for each pair of species in 
the mixture 3 . The basic data are either in tabular form or are given as simple analytical functions of 

temperature or composition. NATA contains 12 methods or options for calculating the cross sections 

— (1.1) — (2,2) 

Trf^y , 7rUy , and B*j (the ratio of cross sections). The options include using the Coulomb cross 

3 Different combinations of the indices l and s are required for higher order terms of the Chapman-Enskog theory. 
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HUtilUi 


3 


section for the electrons and ions plus using exponential potential (ref. 42) and Lennard-Jones (6-12) 
potential (ref. 36) for neutral species in high and low temperature ranges, respectively. The formulas 
used in NATA to compute the transport properties from the collision cross sections are obtained 
from an approximation (ref. 43) to the first-order Chapman-Enskog expressions. These formulas are 
provided in subsequent sections. 


Transport properties of single species . The viscosity p x and frozen 4 thermal conductivity 
Kji of a gas containing a single molecular species are given, to a good approximation, by the 
formulas in chapters VII and VIII of reference 36. First, the viscosity is 


Mi 

= 5 X 10“ 

>6 »sf' 2) 

(17a) 

or 



Mi ~ i 

(2.6693 x 1(T 5 ) 

(17b) 

Further, the viscosity can be written in 

terms of the mean free path \ as 



5tt 

Mi = 

(17c) 

where the gas density is 

p = pi = m l n l 

(17d) 

the mean molecular velocity is 

1 8kT 

U — Uj — \ / 

y irrrii 

(17e) 

and the mean free path is 

V^njTrOf 2) ] 1 x 10 16 


Ai = 

(17f) 


Next, the frozen thermal conductivity can be written in the modified Eucken approximation (ref. 36) 
as 

K fti - K u>i + K in M (18a) 


Here, K^i is the translational component of the thermal conductivity and K\ nty i is the component 
of thermal conductivity resulting from the diffusion of internal excitation energy of the molecules. 
The translational component may be written as 


75 ky/nkT jm % 

tr ’ ,_ 64 ^ 

x 10 16 /J 

(18b) 

= f 1.9891 x 1CT 4 

\ TTs/TjMi 

1 »n< 2 ' 2 > 

(18c) 

15 Mi ^univ 
“ 4 Mi 


(18d) 

_ 3 Mi(^p,i )tr 

” O 


(18e) 


4 See appendix A for various definitions.' 
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or, in terms of the mean free path A^, as 

57T pu\ 


KtT>i 32 Mi 




The internal component of thermal conductivity may be written as 


(18f) 


(18g) 


3 ( C^i 5\ ky/nkT/mi 

8 Uuniv 2 J nU (hD 

6 U 

x 10 16 /J 

(18h) 

(6.3605 x 10~ 5 ) | 

v J Mi 

( Cp.i 3 \ 

^univ 2 / 

(18i) 

Pi ( pDa ^ ( n N 

i (Cp ’ ,)uit 


(18j) 

57T f PiDii\ 

32 Mi V Mi / 


(18k) 


where the internal specific heat has been introduced from the relation 


(^p,i)int — (Ctv)int — 


:Ru 


— n . _ _ y> 

— ^ v.r 2 1 Lu 


(181) 


and the coefficient of self diffusion D u from the first Chapman-Enskog approximation is as follows: 


3 y/irrriikT /1\ 


x 10 16 


= (2.6280 x 10^ 3 ) 


»V t3 /m i 

u r 


(18m) 


(18n) 


The coefficient of self diffusion D l% must be regarded as somewhat artificial. It is more correct to 
regard it as a limiting form of the coefficient of binary diffusion. 

Using equations (18a), (18d), (18j), and (181), the frozen thermal conductivity may now be 
expressed as 


RunivPi 

[15 . 

[t + < 

( P^ii ^ 

f Cp,i 

_5\1 

Mi \ 

^ Mi / 

\ ^univ 

2yJ 


or, from equations (18a), (18e), and (18j), 


£ 

ii 

Ss|5 

2 (^P,i )tr + 

f pDii ^ 
l Mi ) 

|(Cp. 

t)int 

Mi 
” Mi 

2 (^u,i)tr + 1 

f pD n N 
l Pi ) 

|(C., 

,z ) int 


(19a) 


(19b) 

(19c) 
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From equations (17a) and (18m) the factor (pDa/pi) is related to the collision cross sections 7 
through the following relation 


pDu = 6 7 rftg’ 2) 

w 5 7rng ll) 


(20) 


— (2 2 ) — (1 1 ) 

The ratio 7rfi”’ ; /7rf)h in equation (20) is a very slowly varying function of temperature T\ hence, 
pDa/ni is very nearly constant. This factor, appearing in equations (18j), (18k), (19a), (19b), and 
(19c), has a value close to 1.32 for the Lcnnard- Jones potential over a wide temperature range 
(ref. 44). If this factor is approximated by unity, equations (19b) and (19c) for the frozen thermal 
conductivity reduce to the form 



m 

i Mi 

"3 

4“ {Cpj ) jjjt 

(21a) 


Mi , 

Mi 


(21b) 


Ik , 
Mi 1 


(21c) 

or, using equation (17c) to express 

the results in terms of the mean free path, 


K 7. 

5n i 
* = 32 l 

p u \ Wc - + ^ A 

, Mi ) V "’’ ! 4 m "V 

(2 Id) 


This is the form for the thermal conductivity derived originally by Eucken. It differs from the 
modified Eucken approximation (cq. (19)) used in the present work by the ratio 


4 (Ci\i)iut 
15 /? univ 



4 (Cty)int 

1 0 /?univ . 


which is about 6 percent in the case of a diatomic gas near room temperat ure, where the only internal 
energy excited is the rotational energy. 

The modified Eucken approximation (eqs. (19)) neglects the effects of inelastic collisions on the 
thermal conductivity. Such collisions introduce a coupling between the translational and internal 
components of the thermal conductivity; this coupling tends to reduce the total frozen thermal 
conductivity below the value (eqs. (19)) predicted by the modified Eucken approximation. The 
effects have been treated in detail for various polyatomic gases near room temperature. It has been 
found (refs. 45 to 47) that the errors in equations (19) may approach 10 to 20 percent when there is 
a rapid exchange of energy between the internal and translational states through inelastic collisions, 
as is normally the case for rotational excitation in low-temperature polyatomic gases. However, the 
errors become smaller when the exchange is less rapid, and are negligible when 20 or more collisions 
are required for the exchange of energy between internal excitation and translation. 

The effects of inelastic collisions on the thermal conductivity of high-temperature air have 
apparently never been treated in detail. However, it appears that the inelastic cross sections 
should, in general, be small enough to make errors in the modified Eucken approximation (eqs. (19)) 
negligible for air at temperatures greater than about 1000 K. At temperatures below 1000 K, it 
may be desirable for more accurate calculations to correct equations (19) for the effects of inelastic 
rotational-translational energy transfer using the analysis of references 45 to 47. However, this 
correction has not been included in the values presented here. 
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The collision cross sections n or 7 rf are, in general, the weighted averages of the cross 
sections for collisions between species i and j. These collisions have been defined (refs. 7 to 9) as 

= So° So cxp(-7 2 )7 2s '+ 3 (l - cos^ x)47r<Tjj sin y d\ dy 
lJ fo X So cxp(-7 2 )7 2s + 3 (l - cos' x) sin X d\ d-y 

where <7 tJ = <7ij{x,g) is the differential-scattering cross section for the pair (i.j), x > s the 
scattering angle in the center-of-mass system, g is the relative velocity of the colliding particles, 
and 7 = ([m;mj/2(mj + m.j)kT] X/2 ^jg is the reduced velocity. For collisions between the similar 

species, equation (22) yields irfljj required in equations (17) and (18). Various combinations of the 
indices / and s are required for higher order terms of the Chapman-Enskog theory. For the order 
considered here, only 7T<7^ 7rfl|j’ 2 \ and ttH [ 2,2) are needed. 

jq (2 2) 

The collision cross sections 7rfih ! ; and 7rQ- ■ 1 ' employed here are the same as those used by Yos 
(refs. 7 to 9). The cross sections for the neutral species N, N 2 , NO, O, and 0 2 were taken from 
the tabulations of Yun, Weissman, and Mason (ref. 48) for temperatures up to 15 000 K. Above 
15 000 K, the cross sections for atomic N were obtained by extending Yun’s calculations to 30 000 K 
using the same input data and techniques that were used in his work. The cross sections for the 
remaining species N 2 , NO, O, and 0 2 were extrapolated to 30000 K assuming the same temperature 
dependence as calculated for N. 

For the ionized species, the calculations used effective Coulomb cross sections that were chosen 
to make the computed transport properties for a fully ionized gas agree with the correct theoretical 
results (ref. 49) discussed in references 7 and 9. The specific formulas used in the calculations are 


7rn^ ,2) = 1.29Q,. x 10 16 
7r?iJ 2 / 2) = 1.36 Z 4 Q c x 10 16 

= Trfi} 1 /^ = 0.795 Z 4 Q c x 10 16 
where Z = 1 for singly ionized species and 


(for electrons) 
(for ions) 
(for electrons and ions) 


(23a) 

(23b) 

(23c) 


Qr = 


(kT)' 


In A cm 2 


(23d) 


The shielding parameter A is defined as 


A = 


9 (fcT) 3 16(fcT) 2 

47re c rc c + e 4 2 / 3 


1/2 


(23c) 


209 x 10-2 (ipy + 1 52 (dy 


2 / 3 ] ! /2 


(23f) 


where T is the temperature in kelvins and p e = n ( k T is the electron pressure in atmospheres. 6 In 
equation (23d), e = 4.8 x 1CT 10 esu is the electron charge. Equations (23a) to (23f) are applicable 


The collision cross sections designated here as TrClfj*' 1 are the same as given in reference 3G and as 

given in reference 39. 

For flows with thermal nonequilibrium, p t > = n c kT e should he used to obtain the electron pressure. 
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only for electron pressures below the limiting value given by 


Pem = 0.0975 



(23g) 


for which the shielding factor In A is equal to 1 in equation (23d) (ref. 49), and should not be used 
for electron pressures above this limit. The limiting electron pressure, penu from equation (23g) 
is plotted as a function of temperature in figure 5. The electron pressures encountered in typical 
aerospace applications should fall well below this limit. 

In the present report, the transport properties of the charged species (i.e., ions and electrons) are 
provided for the limiting electron pressure p e = p em at which the shielding factor In A = 1. For any 
other electron pressure, it is necessary to correct the tabulated transport properties of the charged 
species according to the formula 


Pi(Pe) 

Pi(pem) 


A^ ? ;(pe) ^tr, i{P?) _ -^int ,i(pfc) _ 1 

^/,i(Pe m) ^t.r,i(Pern) ^int,i(Pem) l n A(Pe) 


In 


2.09 x 10 


-2 


\ 1000p ( 


1/4 


V + 1.52 

) VlOOOpe'V 


8 / 3 ] 


(24a) 


// e 1 ) 

Similarly, the collision cross sections it j (for the pair of species for which both are ions or electrons 
or a combination of the two) for any other electron pressure p e may be obtained from the values 
provided herein for p e = p em by employing the relation 


TTf (Pem) 


A (Pe) 



2.09 x 10~ 2 


, 1000 p^ 4 


4 

+ 1.52 



8 / 3 ] 


(24b) 


Equation (24b) is also applicable for i = j (single species). 

In calculating the contribution of the internal energy states to the thermal conductivity 

for the atomic species N and O, the diffusion cross section in equation (18h) or (18i) has 

been set equal to the corresponding charge-exchange cross section for the atom and atomic ion 
(refs. 8 and 50). As discussed in reference 7, this approximation allows for the effects of excitation 
exchange in reducing the contribution of internal energy states to the thermal conductivity in a gas 
of identical atoms. For consistency, the same approximation has also been employed in the tabulated 


cross sections and self- diffusion coefficients for these species. 

The individual species viscosities and thermal conductivities computed using equations (17) 
and (18) have been curve- fit ted herein as a function of temperature by employing the following 
relations: 

m = [exp(C w )] rK '■> T +^ t ] (25) 




(26) 


The curve- fit coefficients in equations (25) and (26) are given in tables V and VI for the 11 
chemical species. These coefficients yield values of the viscosity and frozen thermal conductivity 
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of ionic species at the limiting electron pressure (eq. (23g)) and should be corrected for any other 
electron pressure (or number density) by using equation (24a). 

Similar to the equilibrium constant calculations from equation (5d), the thermal conductivity 
should be evaluated from the logarithmic form of equation (26). Also, all the five significant digits 
for the coefficients appearing in table VI should be used. 

Figure 6 displays typical results from the viscosity curve fit (eq. (25)) of equations (17) for 
some of the neutral and charged species. The frozen thermal conductivity curve fit (eq. (26)) of 
equations (18) is shown in figure 7 for the same species. Figure 7 and table V show that fourth-order 
curve fits are needed for the thermal conductivity. 

Transport properties of multicomponent mixtures . Rigorous kinetic theory formulas that 
have been derived directly from a solution of the Boltzmann equation using the classical Chapman- 
Enskog procedure (refs. 36 and 51) are available for obtaining the transport properties of a gas 
mixture from the molecular constituent species. In the first Chapman-Enskog approximation, 
formulas for both the viscosity and translational component of thermal conductivity K tr of a gas 
mixture are of the general form 



(27) 


where x i is the mole fraction of the zth species, NS) is the total number of species present in 
the mixture, and the matrix elements Aij can be expressed in the form 


A ij = Aji = -XiXjdij + 6ij | 'xiAi + ^ Xixiau j (28) 

where a tJ = a tI and l) L j is the Kronecker delta. Elements A , and are defined subsequently. The 
superscript 1 on /i or Atr indicates that equation (27) is the first Chapman-Enskog approximation for 
the transport property. Further details for obtaining the transport properties by employing the first 
Chapman-Enskog approximation are given in reference 43. In principle, the problem of calculating 
the transport coefficients for a given mixture consists of two parts: first, the determination of 
the collision cross sections nClij for all possible pairs of species and second, the evaluation 

of the Chapman-Enskog formulas. The amount of computation required to evaluate the mixture 
transport properties is greatly reduced if approximations to the complete Chapman-Enskog formulas 
are employed. Reference 7 has provided approximate formulas for the transport properties based 
on the relations developed in references 52 to 54. References 13, 14, and 55, and more recently, 
references 15 and 56 have also provided approximations to the Chapman-Enskog formulas. These are 
apparently the most satisfactory of the many simplifying approximations for the mixture viscosity 
and thermal conductivity that have been suggested by various authors. However, effects of the 
elements a^j in the Chapman-Enskog formula are completely neglected in Brokaw’s approximation 
(ref. 53), so that this approximation always gives too large a value for the transport properties. In 
the Buddenberg- Wilke (refs. 13 and 55) and Mason-Saxena (ref. 14) formulas, the effects of these 
elements are accounted for by means of a single empirical constant, which is assumed to be the same 
for all gas mixtures. The approximation used by Peng and Pindroh (ref. 52) represents an attempt 
to take account of the nondiagonal elements explicitly to the first order at the expense of a somewhat 
increased calculational effort. Armaly and Sutton (refs. 15 and 56) neglected the nondiagonal matrix 
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elements A tJ 
defined as 


in a manner similar to Brokaw (ref. 53). However, they did not force the value of A*^ 


n (2 - 2) 

A*. - 

lj n (u) 

ij 


(29) 


to be equal to 5/3 and 5/2 in their approximations for viscosity and thermal conductivity, 
respectively. 7 They assigned different values to A*j for ion-atom and neutral atom-molecule 
interactions. From the computer time, storage, and simplicity point of view, references 13 and 14 
appear to be adequate for nonionized gas mixtures, whereas references 15 and 56 are useful for 
computing the viscosity anti translational component of thermal conductivity for an ionized gas 
mixture. 

In all the approximations to the Chapman-Enskog formulas for viscosity and translational thermal 
conductivity discussed thus far, the transfer of momentum or energy from one species to another 
by collisions has been either neglected or has been accounted for by an empirical constant. This 
transfer process, which is represented by the nondiagonal elements a t j in the Chapman-Enskog 
formula (cq. (28)), has the effect of making the less conductive species in the mixture carry a larger 
fraction of the transport. This process, therefore, reduces the overall conductivity of the mixture 
below the value it would have if the transfer process were neglected. 

In reference 43, Yos obtained approximations to the Chapman-Enskog formulas that account for 
the effects of the aforementioned transfer process between different species. These approximations 
reproduce the results of the first Chapman-Enskog formula (eqs. (27) and (28)) to within a fraction 
of a percent for air at 6000 K and for other cases considered in reference 43 and are simpler to use 
than the latter. Based on the relations developed by Yos (ref. 43), the following formulas may be 
used to compute the mixture viscosity and translational component of thermal conductivity: 
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Here, NS is the number of species in the gas, x t is the mole fraction of the ith species, and a a v is an 
average value of the nondiagonal matrix elements contained in equation (28) and is defined as 
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For the viscosity, the quantities a,j and B tJ in equations (31a) and (31b) are defined as 
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7 With A\- = 5/3 and 5/2, A tJ becomes identical to zero in the Chapman-Enskog formula. 
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For translational thermal conductivity, the above quantities are defined as 
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In these equations, Mi is the molecular weight of the ith species, N a is the Avogadro number, and 
k is the Boltzmann constant. The remaining quantities are defined 
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and the collision integrals n are weighted averages of the cross sections defined in equation (22). 

The mixture frozen thermal conductivity K j employed in defining the frozen Prandtl and Lewis 
numbers 8 can be obtained from the modified Eucken approximation (ref. 54) 


Kj — K tr + ^int (37) 

where Ki r is the translational component of the thermal conductivity given b} r equations (30) to (3G) 
and K mi is the component of thermal conductivity resulting from the internal excitation energy of 
the molecules, (ref. 7) given by 
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If the approximations of references 52 to 54 are employed in place of the more exact formulas 
given by equations (30) to (36), the following approximate formulas for the mixture viscosity and 

8 See appendix A for various definitions. 
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thermal conductivity may be used (ref. 7). For the mixture viscosity 
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For the translational component of thermal conductivity in a mixture, 
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is given by equation (35) and is defined as 
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The approximations of equations (39) and (40) may be valid only for the range of conditions for 
which they have been developed and not for general application because of the approximate analysis 
of the nondiagonal matrix elements Ajj of equation (27). Generally, these approximate formulas 
give good results for nonionized or weakly ionized flows with a savings of about a factor of two in 
computation time compared with the more accurate formulation of equations (30) to (33), since it 
is no longer necessary to evaluate the nondiagonal matrix element a av in equation (31a). 

A further savings of about another factor of two in the transport property calculations may be 
obtained if a in equation (40b) and the factor 
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equations (38) to (40) depend on the same parameter i n ^e denominator. With these 
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approximations, equations (38) to (40) may be written in the form 
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where the partial gas density of the ith component is 


Pi = mini 

the mean molecular velocity of the ith component is 

U{ — \J SkT/ixmi 

and the modified mean-free path of the ith component in the mixture is 
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In the above equations, n is the total particle density in particles/cm 3 and = (Cp Onit- The 

modified mean- free path X x given by equation (41f) is not the mean distance between molecular 
collisions. It differs from the actual mean-free path A* of the ith component in a mixture defined in 
reference 6 as 
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by the factor 2(1 + in each term of the summation. The modified mean-free path \ x 

represents the mean distance travelled by species i before losing its average excess momentum 
relative to the surrounding gas medium through head-on collisions with species j. (See Hansen, 
ref. 57.) From equations (37), (41b), and (41c), the expression for the frozen thermal conductivity, 
Kj, of the mixture can be written as 
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in close analogy to the Eucken expression (eq. (21d)) for a single-component gas. Equations (41a) 
to (41i) represent the predictions of the elementary kinetic theory for the viscosity and thermal 
conductivity of a gas mixture. (See refs. 36 and 58.) Although less accurate than the approximation 
given by equations (37) to (40), they may still be adequate for many applications, particularly in 
view of the current uncertainty in our knowledge of the collision cross sections for high-temperature 
gases. 

The binary diffusion coefficient needed to obtain the binary and multicomponent Lewis 
numbers (refs. 17, 18, and 36) is obtained from the complete first Chapman-Enskog approximation 
(ref. 36): 
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Here, p is the pressure in atmospheres. Equation (42b) for D(j has been curve fitted in the present 
work by using the expression 
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where T is the temperature in kelvins. The curve-fit coefficients of equation (42c) arc given in 
table VII for the different interactions in the 11-species air model. 

The values of D tJ obtained from equation (42c) with the curve-fit coefficients given in table VII 
are for the limiting electron pressure p e m (cq. (23g)). 9 If the pair of interacting species are both ions, 
both electrons, or any combination of the two, then D tJ must be corrected for the given electron 
pressure as follows: 


Di i {pf ) — Djj(p c in) 


In 


2.09 x 10“ 2 


(■ 


V iooop, 


p-y 


+ 1.52 


, \ 8/3' 

V 1000p ( 1/4 ) 


( 


(42d) 


To obtain the viscosity p and the frozen component of thermal conductivity Kj of a gas mixture 
either from the more exact equations (30) to (38) or from the approximate relations in equations (39) 
to (41), the binary collision cross sections TrfiJ- and their ratios are needed in equations (34) to (36). 
These cross sections, defined by equation (22), arc the same as those used by Yos (refs. 7 to 9). In 
addition to the formulas given by equations (23) for the Coulomb collision integrals, the following 
relations are employed: 


7rfi|. 1 / ’ 1) = 0.795 Z 2 Q c x 10 1G (43a) 

*ng' 2) = 1.29 Z 2 Q c X 10 16 (43b) 

where the subscripts e and / represent electrons and ions, respectively, and Z = 1 for singly ionized 
species, and Q c is defined in equation (23d). 

The NATA (Nonequilibrium Arc Tunnel Analysis) code (refs. 39 to 41) employed to obtain the 

collision cross sections 7 rfl^ ‘^contains default provisions for estimating some cross sections if they 
are not specified explicitly in the built-in data base or the input. The defaults are summarized as 
follows: 

1. If both species are ions, the Coulomb cross sections given by equations (23) and (43) arc 
used. 

2. If one species is neutral and the other is ionized, the formula 

7T = A {Ls) T~ 0¥ (44a) 

is employed with the constants A defined in the code. 

3. If both species are neutral and are not alike, the cross sections are estimated using the 

simple mixing rule 
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9 For flows with thermal equilibrium, p ( = n c kT may be used to obtain the electron pressure. For thermal nonequilibrium 
conditions, the electron temperature T c should be employed instead of T in the previous equation of state. 
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The built-in data in NATA specify steps for calculating the cross sections for the like-like 
interactions of ten species (e”, N 2 , O 2 , N, 0 5 NO, N0 + , N + , 0 + , and Nj) and for those unlike 
interactions for which experimental or theoretical cross sections arc available in the literature. The 

cross sections for Oj are the same as those for Nj in these calculations. NATA contains 12 

n n (2 2) 

methods or options for calculating \ 7rfl ■' , and B*-; these methods are described in detail 

in references 39 to 41. Reference 59 provides a comparison between the theoretical calculations 
based on these cross sections and the experimental values for the thermal conductivity of nitrogen 
up to temperatures of 14000 K; fairly good agreement between the two is shown. The accuracy 
of the calculated transport properties is largely determined by the accuracy of the input data for 
the cross-section integrals. Reference 7 contains a detailed discussion of the values employed in the 
present study, and a brief review of these values is provided in the following paragraph. 

Similar to the collision cross sections for single-species transport properties described previously, 
the values of 7rfi} ■’* Tor mixture transport properties for the atomic and molecular interactions have 
been obtained from the calculations of Mason (ref. 48) for temperatures up to 15000 K. The cross 
sections for the interactions not tabulated by Mason were assumed to be approximately the same 
as the other neutral-neutral interactions. For temperatures above 15000 K, the cross sections for 
the atomic and other neutral-neutral interactions were obtained by extending Mason’s calculations 
up to 30000 K by using the same potential functions that were used in his work (ref. 60). Since 
neutral species are not important at high temperatures, the effect of this procedure on the overall 
transport properties should be negligible. Further, there appear to be insufficient data available on 
the ion-neutral interactions. These interactions have been calculated from the approximate potentials 
given by Peng and Pindroh (ref. 52). Those potentials were obtained by drawing a smooth curve 
joining the polarization potential at large intermolccular distances to a Morse potential obtained 
from spectroscopic data at short distances. Further details for obtaining the interactions between 
a neutral species and its own ion with resonant and nonresonant ion-neutral collisions are given 
in references 7 and 8. There has been a lot of experimental work on the electron-molecule cross 
sections, whereas the cross sections for electron-atom scattering are somewhat uncertain, as detailed 
in reference 7. For Coulomb collisions between the charged particles, effective collision integrals 
were chosen in reference 7 so as to make the calculated values of the electrical conductivity and the 
electronic contribution to the thermal conductivity agree as closely as possible with the results for 
a completely ionized gas. 

Figures 8 and 9 show viscosity and frozen thermal conductivity values, respectively, of equilibrium 
air at 1 atm obtained by employing the collision cross sections computed here for the constituent 
species. The present calculations employ equation (39) for the mixture viscosity and equations (37), 
(38), and (40) for the mixture frozen thermal conductivity. Also, the collision cross sections used 
are nearly identical to those obtained in reference 9. The mixing laws employed for the viscosity and 
thermal conductivity shown in figures 8 and 9 are accurate for subionization temperatures (less than 
9000 K at 1 atm) only. For ionized flows, more accurate mixing laws of the type given in equation (30) 
should be employed. Figures 8 and 9 also include viscosity and thermal conductivity values from 
other sources for comparison. In these figures, the predictions of Peng and Pindroh (ref. 52) and Esch 
et al. (ref. 61) are based on the Buddenberg- Wilke (ref. 55) type of mixture law and have the same 
level of approximation as the present calculations. The mixture laws employed by Hansen (ref. 37) 
and Svehla (ref. 44) contain a somewhat lower level of approximation. In Hansen’s and Svehla’s 
work, the viscosity is computed by using the simple summation formula for a mixture of hard-sphere 
molecules, whereas a linearized expression with Eucken’s assumption (ref. 36) is used for the frozen 
thermal conductivity. The viscosity values obtained in the present work arc in good agreement 
(fig. 8) with the values obtained by Peng and Pindroh (ref. 52) and Esch et al. (ref. 61), presumably 
because of the similar mixing laws. Hansen’s (ref. 37) predictions, which are based on the Morse 
potential function, are lower, especially at higher temperatures. Svehla’s (ref. 44) predictions, which 
are based on the Lennard-Joncs potential, are also lower than the present results for temperatures 
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higher than 2500 K. Even with the beginning of dissociation of molecular nitrogen at about 4000 K 
(when the dissociation of molecular oxygen is almost complete), Hansen’s values of viscosity are 
not much different from Sutherland’s law. For the frozen thermal conductivity (fig. 9), the present 
values are in agreement with those obtained by Peng and Pindroh (ref. 52) up to temperatures of 
about 9000 K. The values obtained by Esch et ah (ref. 61) deviate from the present values beyond 
6000 K. This deviation may be due to the constant cross-section values that are used for the ionized 
species in the 8000 to 15000 K temperature range in reference 61. Hansen’s predictions of thermal 
conductivity are lower than the other data and are closer to the Sutherland values up to temperatures 
of about 4000 K. Again, the differences between the present computations and those of Hansen are 
presumably due to the somewhat more rigorous mixing laws employed herein. Svehla (ref. 44) has 
not provided values of the frozen thermal conductivity for equilibrium air and, therefore, his values 
are not shown in figure 9. 

There are 121 possible binary interactions for the dissociating air with 11 species. Therefore, 

— (1 l ) ' —(2 2 ) 

121 values of each of the collision cross sections 7r , 7 rfi)- 1 , and the collision cross-section 
ratio Z?* are required to evaluate the transport properties. If the symmetrical equality is used, 
i.e., (i,j) = (j,z), only 66 values of each cross section and cross-section ratio are required. These 
values have been curve-fit in the present study as a function of temperature for the limiting electron 
pressure p em (eq. (23g)) using the following relations 
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Curve-fit coefficients in equations (45) to (47) are given in tables VIII to X for all interactions in the 
11-species air model. Similar to the calculations for the equilibrium constant and thermal conductiv- 
ity, the collision cross sections should be evaluated from the logarithmic forms of equations (45) and 
(46) by keeping all four significant digits for the coefficients given in tables VIII to X. For electron 
pressures different from p em , the formula given by equation (24b) is used to correct the cross sections 
for the ionic species. No such correction is required for the cross-section-ratio parameter B*y 

Figure 10 illustrates some typical curve-fits that were obtained by using equation (42c) and 
equations (45) to (47) with the associated constants. The figure compares the computed values of 
binary diffusion coefficient, collision integrals, and collision- integral ratio with the resulting curve-fit 
for different interaction pairs of neutral and ionized species, including electrons. The collision-integral 
ratio, B *. , is almost constant with temperature as shown in the figure and was fitted with the lower 

— z l 

order curve-fit when possible. The collision integrals * for the charged species pairs show a 
simple T -2 dependence that requires only two curve-fit coefficients, as seen in equations (23). 


Thermal Nonequilibrium 

The transport properties of a multitemperature gas mixture may be obtained by following the 
approaches of references 27 and 62. These references have used equations (39) and (40) for the 
calculation of mixture viscosity and translational thermal conductivity. The collision integrals 
for heavy particles in these equations are evaluated by using the heavy-particle translational 
temperature, whereas those for electrons with any other partner are obtained by using the electron 
temperature. With these modifications, equations (39) and (40a) can be written for a gas mixture 
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where a ZJ is still obtained from equation (40b). The above definition of K* r does not include 
contributions due to electron-heavy-particle and electron-electron collisions. These contributions, 
defined by K e , are given subsequently. The approach outlined here may also be used to obtain the 
mixture viscosity and thermal conductivity from the more exact expressions in equations (30) to (36). 
Equations (38) for the internal thermal conductivity need to be modified for the multitemperature 
formulation. The contributions resulting from the excitation of different internal energy modes 
cannot be lumped together into a single term K mt for such a formulation. Further, equation (37) 
can no longer be used to obtain a frozen thermal conductivity for the mixture. The relation given 
in equation (37) may be used to obtain a frozen thermal conductivity with only the rotational 
mode contributing to the internal energy at the translational temperature. In general, there are 
four components of the internal thermal conductivity, similar to the molecular specific heat (fig. 4). 
Using these components, the kth component of the overall heat-flux vector can be expressed as 
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where K* r is the translational thermal conductivity defined previously by equation (49) and A" ro t, 
K v fo, and K e \ are the rotational, vibrational, and electronic thermal conductivities, respectively, 
associated with these internal energy modes (ref. 63). Also, K e is the thermal conductivity of 
electrons, T is the translational-rotational temperature, T v ^ is the vibrational temperature, T e i is 
the electronic excitation temperature, T e is the electron temperature, x k is the kth component in a 
general orthogonal coordinate system, and the last term is the diffusive component of the heat-flux 
vector. In the diffusive heat-flux component, p 2 *, h t , and V k are the density, enthalpy, and diffusion 
velocity of species i, respectively. Further details of the heat-flux vector q k and other definitions are 
given in appendix A. For a two- temperature model, equation (50a) may be written (refs. 25 and 27) 
as 
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Different components of the internal thermal conductivity in equations (50) can be evaluated 
from equation (38a) or (38b) by appropriate modifications for these components. For example, A" ro t 
can be obtained from 
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for partial excitation of the rotational internal energy mode if the temperature is less than that 
needed to excite the vibrational energy mode (fig. 3). Values of specific heat at constant pressure 
i in equation (51a) can be obtained from the curve-fit relation of equation (13) by employing the 
rotational temperature T ru t if different from the translational temperature T. In equation (51b), 
(C Pli ) ro t is the rotational component of the total C p j. When the vibrational mode begins to excite, 
the rotational mode is fully excited, and equation (51a) becomes 
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Similarly, expressions for K v \\ } with partial and full excitations of the vibrational energy mode may 
be written as 
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Once again, the value of C p j in equation (52a) can be obtained from the curve-fit relation of 

equation (13) by employing T v ^; also, in equation (52b), (C'p.Ovib is the vibrational component 

of the total C p j . The vibrational energy mode is fully excited when the electronic contribution 

becomes significant (as shown in fig. 4 for the specific heat at constant pressure). 

The rotational and vibrational energy modes are almost fully excited at their respective 

characteristic temperatures. Therefore, for rotational and vibrational temperatures greater than 

these characteristic temperatures, equations (51b) and (52b) can be used to obtain K ro t and /\ v jb, 
respectively. - 
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Expressions for K v \ with the excitation of electronic energy mode can be written as 
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In equation (53a), higher electronic degeneracy levels begin to contribute to K c \ (and C p j. etc.) as 
the characteristic temperatures for the electronic excitation of those degeneracy levels are reached. 
The term (C pp ) ( .\ in equation (53b) is the electronic contribution to the total C jKl . The species that 
is left out in the outer summation in equations (53a) and (53b) is the electron. 

Finally, the thermal conductivity for free electrons K c in equations (50) may be obtained from 
the modified form of equation (40a) as follows: 
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This thermal conductivity results from the collisions between electrons and other species, including 
other electrons. 

In equations (51) to (53) for the partial conductivities in thermal nonequilibrium, the diffusion 
rates of the excited species are assumed to depend on the translational (or kinetic) temperature of the 
molecules rather than on the excitation temperature, because the collision frequency of the species 
is likely to depend strongly on the translational temperature. Therefore, the cross sections 
in the denominators of these equations are evaluated at the translational temperature T. Perhaps 
employing some averaged temperature, which would cover the effects of translational temperature 
and excitation on the cross sections, would be more desirable. 


Concluding Remarks 

The present work provides a review of the reaction-rate coefficients and thermodynamic and 
transport properties for an 11-species air model that are needed in analyzing the high-energy flow 
environment of currently proposed and future hypersonic vehicles. The properties not available in 
the literature have been provided for this set of species, and curve fits are given for the various species 
for their efficient computation in flow-field codes. Approximate and more exact formulas have been 
provided for computing the properties of partially ionized air mixtures in chemical and thermal 
nonequilibrium around such vehicles. Limitations of the approximate mixing laws for a mixture 
of ionized species are pointed out, and an electron number-density correction for the transport 
properties of the charged species is given. This correction has generally been ignored in the aerospace 
literature. 

The work presented here uses the best estimates of available data needed to compute properties of 
the 11-specics air model. However, there is need for improved data, especially for air in thermochem- 
ical nonequilibrium. There is a considerable degree of uncertainty about reaction-rate coefficients 
at high temperatures. For the multitemperature kinetic models, the theoretical basis needs to be 
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developed substantially and the experimental data base needs to be expanded considerably. The 
virial-coefficicnt and partition-function approaches for obtaining the thermodynamic properties are 
equally accurate at high temperatures if the effects of nonrigidity of the rotor and anharmonicity 
of the oscillator are included in the two methods. The curve fits provided here are for the thermo- 
dynamic properties based on the virial-coefficient formulation with these high-temperature correc- 
tions. These values seem adequate but need further verification. For transport properties, the input 
data for obtaining the collision cross sections at high temperatures require further study. 



Appendix A 

Heat Flux, Frozen and Total Prandtl and Lewis Numbers, and Associated Definitions 

The fcth component of the overall heat-flux vector q k for the dissociating and ionizing air model in a 
multitemperature formulation can be expressed as follows: 10 
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where the various symbols are explained in the main text after equation (50a). For thermal equilibrium 
conditions, equation (Al) may be written as 
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Also, under thermal equilibrium conditions, 


Ktr = K* + K t 
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where Ki r is given by equation (40a), K* T is obtained from equation (49), and K e is evaluated from 
equation (54). The last terms in equations (Al), (A2), and (A3) represent the diffusion contribution to 
the heat-flux vector. The diffusion mass flux of species i, J k , is related to the diffusion velocity V k through 
the relation 

Jf = Pi v k (A7) 

or, in terms of the concentration gradients, may be expressed (refs. 64 and 65) as 


j? = - M 


NprJ 


‘ NS \ 

r 9C i V- a » 9 C t 




l dx k "p u dx k 
¥*' 


(A8) 


where pressure and thermal diffusion have been neglected and where 


NS C . /NS c 


^-ES/E- 


jTi M i / pi M j Le f,ij 


(A9) 


A 6,7 = Le, — 


Mi 

M 


/ M \ NS 


(A10) 


10 The heat-flux vector is positive for heat flow in the positive coordinate direction. 
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The multicomponent Lewis number Ljjj and the binary Lewis number Le j ^ are defined, respectively, as 

L/ij = pCpjDij/K j (All) 


and 


L e fSj pCp j D t j /Ay 


(A12) 


where D tJ is the multicomponent diffusion coefficient that can be obtained from the binary diffusion coefficient 
following the approach of reference 65. 

For binary diffusion, equation (A8) may be simplified to (ref. 64) 




(A 13) 


If the species j is assumed to be some effective mean species for the gas, then using equations (A7) and (A13) 
in equation (A3) gives 

(AH) 


<1 


1=1 


or, for chemical equilibrium flow applications, 


d T 


NS 


dC; \ OT 


q ~ ~ Kf dx k ~ I Dijhi lrr TuT 


i = 1 


OT J dx k 


or 


H dx k 

where K is the total effective thermal conductivity defined as 


NS 


„ „ L 0C t 

T — Ay + P ^ Djjhi ^ 

i— 1 


or 


K — Kj + K r 

with the reactive contribution to the thermal conductivity K r defined as 


« r = pE D ‘i l ‘igf 

i= 1 


(A15) 
(A 16) 

(A 17) 
(A18) 

(A19) 


for chemical equilibrium conditions. 

An alternate definition of the reactive thermal conductivity 11 K r for a flow in chemical equilibrium is given 
in reference 7 as 


NIR 

1=1 


(Ahi/ R au \ v T) 2 


NS NS 

T, - a i,l)/ x i] T [($,* _ a i,l) x i - (0j,l - a j,l) x i] A 
i= 1 j = 1 


(A20) 


( 1 ) 

V 


11 The formula for K r given here already includes the effects of ambipolar diffusion (ref. 66} on the reaction conductivity of an ionized 
gas. . ... . 
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where k is the Boltzmann constant, NIR is the total number of independent reactions in the system, NS is 

NS 

the total number of species in the system, A h[ — Y {@il ~ a il)hi is the heat of reaction per g-mole for 

i = 1 ' ’’ 

the Ith reaction, xi is the mole fraction of species z, A.j^ is defined by equation (34), and otjj and are 
the stoichiometric coefficients for reactants and products in the reaction given by equation (1). The thermal 
conductivities K r given by equations (A 19) and (A20) are equivalent under conditions of chemical equilibrium. 
Similar to relations (A17) and (A18), the total specific heat at constant pressure may be defined as 


or, from equation (9), 


or 


with 


ft. 

II 

(A21) 


(A22) 

Cp = Cpj + C Pr 

(A23) 

C -(S '••§), 

(A24) 


where C Pr is the total contribution to the specific heat at constant pressure. 
The specific heat at constant volume C v may be obtained from 


C v 

where 

NS 

e = E 

?'=i 

This equation may be used with equation (A21) to 

7 

Using the various definitions for the thermal conductivity and the specific heats at constant pressure, the 
frozen and total Prandtl and Lewis numbers usually employed in dimensionless heat- flux calculations may now 
be defined. 12 

Generally, the frozen values of specific heat at constant pressure and thermal conductivities arc employed 
in flow-field calculations. Associated heat flux from these calculations is, accordingly, expressed in terms of the 
frozen Prandtl and Lewis numbers. Alternatively, one may also use total values of specific heat at constant 
pressure and thermal conductivity in flow-field computations. In this case, the heat flux may be expressed in 
terms of the total Prandtl and Lewis numbers. However, the total values of C p and K , as well as those of the 
total Prandtl and Lewis numbers, can be used only with calculations involving thermal equilibrium. This is 
obvious from the definitions of K and C p given by equations (A18) and (A23), respectively. 

The frozen Prandtl number and Lewis numbers are defined as 


■(ft). 


(A25) 


(hi Vi! Pi] 


(A26) 


obtain the ratio of specific heats 7 from 

= Cp/Cv 


(A27) 


N PrJ = CpjufKj (A28) 

Le f,ij = pCp f D ij/ K f (A29) 

12 Kenneth Sutton of NASA Langley Research Center provided some of these definitions in an unpublished memo. 
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where C P f and Kj are defined by equations (9) and (37), respectively. Generally, the subscript / is not used 
in the literature to denote the frozen Prandtl and Lewis numbers as is done here. 

The total Prandtl and Lewis numbers (consisting of both frozen and reactive components) are defined by 
employing the total values of C p and K as 

N Pr = C p p/K (A30) 

Le lJ - pC p D i: j/K (A31) 

where C p and K are defined through equations (A22) and (A18), respectively. Reference 37 gives similar 
definitions of the total Prandtl and Lewis numbers. The frozen values defined by equations (A28) and (A29) 
are denoted as partial values in reference 37 and are denoted as Pr 1 and Le f therein. 


| 

| 
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Appendix B 

Procedure For Obtaining the Heats of Formation at 298.15 K 

As mentioned in the main text, various sources provide tabulated values of the species specific enthalpy and 
free energy based on a reference temperature of 0 K. To transform these thermodynamic values to the 298.15 K 
reference, the heat of reaction is required at the new reference temperature. The general method employed for 
this transformation involves calculation of the heat of formation for the reaction 

Reactants Products 

53 ot it rXi — ♦ X) ( B1 ) 

t= 1 i= 1 

For illustration, the following reaction for the formation of atomic nitrogen, N, may be considered: 

^N 2 -» N (B2) 

The heat of formation of a substance is defined as the change in enthalpy created by its production. Since 
the enthalpy is a point function of temperature and is independent of the path used to arrive at that point, 
the reactants may be cooled from the reference temperature of 298.15 K to 0 K, thus allowing the formation 
reaction to occur at the lower temperature. The products may then be heated to the reference temperature of 
298.15 K. Through this process we can obtain the heat of formation of atomic nitrogen, N, at 298.15 K by the 
following relation: 

(A^n)t=298.15 = (A^n)f= 0 + A/li + A^2 (B3) 

where 

Reactants 

A/ii = 5Z Qt ’ r [Oi)r=o - Oi)r=298.i5] (B4a) 

i = 1 

= 2 [( /i N 2 )t=0 - (^N 2 )t= 298.15] (B4b) 

Products 

A/i 2 = 53 @h r 298.15 “ Oz)r=o] (B5a) 

i=l 

= [(^n)t=298.15 - (^n)t=o] (B5b) 

The reference species used here as reactants in the formation reactions are N 2 , O 2 , and e - . Care should be 
taken in selecting the reaction for the formation of a particular compound. The reactions used for the formation 
of each substance, along with the heats of formation at 0 K and 298.15 K, are listed in table Bl. The heats 
of formation at 0 K are taken from reference 38. The procedure for obtaining the heat of formation for any 
substance at 298.15 K from the corresponding value at 0 K is illustrated in figure Bl. 
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Table Bl. Reactions for the Formation of 


Substance 

Reaction 

N 

1 

2 

n 2 

-» N 

O 

1 

2 

0 2 

— 0 

N+ 

1 

2 

n 2 

-> N + + e~ 

0 + 

1 

2 

02 

0 + + c~ 

NO 

1 

2 

N 2 + \ 0 2 NO 

N0+ 

1 

2 

N 2 + ^ 0 2 -» N0++ e" 

N 2 

n 2 -* 

N+ + e- 

01 

O2 ~ > 

0 2 + + e- 


Substances With Heats of Formation 


(Ahj ) T = 298.15^ 
kcal/g-molc 

(A/if ) r=0 , 
kcal/g-mole 

112.973 

112.529 ± 0.024 

59.553 

58.984 ± 0.024 

449.840 

447.694 ± 0.1 

374.949 

373.024 ± 0.024 

21.580 

21.457 ± 0.04 

236.660 

235.180 ± 0.2 

360.779 

359.298 ± 0.01 

279.849 

278.370 ± 0.2 




Reactants 
at 298.15 K 


Reactants 

Ahi = a j r [(hj) T=0 


( h i)T=298.15] 


Reactants 
at 0 K 


(Ah j f ) T=298.15 

is unknown 


(Ah j ) o 

is known 
from 

tabulated 

data 


T 


T 


Products 
at 298.15 K 


Products 


Ah 2 = £ P j, r [(hj) t=298.15 ^Pt=(J 


i = 1 


Products 
at 0 K 


Figure Bl. Flow chart for obtaining heat of formation at 298.15 K from corresponding value at 0 K. 
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Appendix C 

Sample Program To Evaluate Thermodynamic Properties From Polynomial Curve 
Fits 

This appendix is a sample Fortran subroutine which evaluates the species specific heats and enthalpies for an 
11-species air model using polynomial curve fits as functions of temperature. Five temperature ranges are used 
for each species for temperatures between 300 K and 30 000 K. Properties evaluated near the temperature range 
boundaries are smoothed by linearly averaging the polynomial coefficients to assure continuous derivatives. The 
subroutine may be easily modified for different needs. 

SUBROUTINE THERMO (T, CPI , HI) 

C Computes enthalpy and specific heat for il species by 
C approximating polynomials. Polynomial coefficients are stored in 
C arrays A1 to A6 and are linearly averaged at the temperature range 
C boundaries . 

C 

C input: T temperature, K 

C output: CPI specific heats of the species, cal/g-mole-K 
C HI enthalpies of the species, cal/g-mole 

DIMENSION Al(ll,5) ,A2( 11,5) ,A3( 11,5) ,A4( 11,5) ,A5( 11,5) ,A6 (11 ,5) 

DIMENSION P(6) ,C0EF(11 ,5,6) 

DIMENSION CPI (11) , HI (11) 

EQUIVALENCE (A1,C0EF) 

C Universal gas constant, cal/g-mole-K 
DATA UNIR /l. 987/ 

C 

C Coefficients are input for five temperature ranges 
C 

K=4 

L=5 

IF(T.GT. 15500.) GO TO 20 

K=3 

L=4 

IF (T . GT . 6500) GO TO 30 

K=2 

L=3 

IF(T.GT. 1200 . )G0 TO 40 

K=1 

L=2 

PA=1 . 0 
PB=0 . 0 

IF(T . LE. 800 . )G0 TO 50 
PB= ( 1 . / 400 . ) * (T-800 . ) 

PA=1 . 0-PB 
GO TO 50 
40 CONTINUE 
PA=1 . 0 
PB=0 . 0 

IF(T . LE. 5500 . )G0 TO 50 
PB=(1 . / 1000 . ) * (T-5500 . ) 

PA=1 . 0-PB 
GO TO 50 
20 CONTINUE 
PA=0 . 0 
PB=1 . 0 
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IF (T.GE. 25500. )G0 TO 50 
PA=1 . 0 
PB=0 . 0 

IF(T. LE. 24500. )G0 TO 50 
PB=0 . 001* (T-24500 . ) 

PA=1 . O-PB 
GO TO 50 
30 CONTINUE 
PA=1 . 0 
PB=0 . 0 

IF (T . LE . 14500 . ) GO TO 50 
PB=0.001*(T-14500.) 

PA=1 .O-PB 
50 CONTINUE 
C 

T2=T*T 
T3=T2*T 
T4=T3*T 
T0V=1 . 0/T 

C Compute properties for 11 species 
DO 65 1=1,11 
DO 60 J=1 ,6 

60 P ( J) =PA*COEF (I ,K, J) +PB*COEF ( I , L , J) 

HI (I) =UNIR*T* (P C 1) +0 . 5*P (2) *T+P (3) *T2/3 . +0 . 25*P (4) *T3 
1 +0 . 2*P (5) *T4+P (6) *TOV) 

CPI (I) =UNIR* (P ( 1) +P (2) *T+P (3) *T2+P (4) *T3+P (5) *T4) 

65 CONTINUE 
RETURN 
END 
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Table I. Third Body Efficiencies Relative to Argon 


Efficiencies relative to argon of 


j— NS),i 

l>i 

2.1 

3.1 

4.1 

5.1 


n 2 

(i = 2) 

0 

(i = 3) 

N 

(i = 4) 

NO 
(i = 5) 

2 

25 

1 

1 

2.5 

1 

0 

1 

1 ' 

20 

20 

20 

1 

0 

0 

0 

0 

0 

0 

0 


NO+ 
(i = 6) 
0 
0 
0 
0 
1 


(i = 7) ( i 

0 

0 

0 

0 

1 


(< = 9) (t = 10) 


*Extension for the 11 species and 5 catalytic bodies is based on the work of reference 17. 
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Reactions and reaction rates for numbers 1 to 7 are from Blottner’s (ref. 17) 7-species chemical model and those for numbers 8 to 20 are from 
reference 21. 


Table III. Curve-fit Coefficients for Equilibrium Constant K eq at Different Number Densities 
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Table IV. Constants for Polynomial Curve Fits of Thermodynamic Properties for a Reference State* of 298.15 K 
and 1 atm Pressure With N 2 , 0 2 , and e“ as Reference Elements (300 K < T < 30000 K)t 
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with T re f _ 298 K, C p q — 5.44 cal/g-mole-K, C p Q 2 — 7.00 cal/g-mole-K, C p >[ = 4.97 cal/g-mole-K, and C p ^ 2 = 6.96 cal/g-mole- 



Tabic V. Constants for Viscosity* Curve Fits 
(1000 K < T < 30000 K)t 



* Viscosity is obtained in g/cm-sec. 

f For temperatures less than 1000 K, Sutherland’s viscosity law for air may 
be used for all species. 

t The charged species viscosities given here are for the limiting electron 
pressure p e m (cq. (23g)). For different electron pressures, these values should 
be corrected by the formula given in equation (24a) in the main text. 




Table VI. Constants for Frozen Thermal Conductivity* Curve Fits 
(1000 K < T < 30 000 K)t 


Species*!- 

A Kf, 

B *u 

° K f,i 

Dk i, 

E «f.i 

n 2 

0.03607 

-1.07503 

11.95029 

-57.90063 

93.21782 

02 

.07987 

-2.58428 

31.25959 

-166.76267 

321.69820 

N 

0 

0 

.01619 

.55022 

-12.92190 

O 

0 

0 

.03310 

.22834 

-11.58116 

NO 

.02792 

-.87133 

10.17967 

-52.03466 

88.67060 

NO+ 

-.06836 

2.57829 

-35.72737 

219.09215 

-519.00261 

e~ 

0 

0 

.00032 

2.49375 

-27.89805 

N+ 

0 

0 

.03088 

2.06339 

-31.51368 

0+ 

-.04013 

1.32468 

-16.22091 

89.96782 

-208.57442 

N 2 + 

0 

-.03723 

.84192 

-3.59040 

-18.65620 

0 2 + 

-.08373 

2.75459 

-33.74529 

185.13274 

-401.50753 


* Thermal conductivity is obtained in cal/cm-sec-K. 

f For temperatures lower than 1000 K, Sutherland’s law for thermal conductivity of air may 
be used for each species. 

J The charged species frozen thermal conductivities given here are for the limiting electron 
pressure p e m (eq. (23g)). For different electron pressures, these values should be corrected by 
the formula given in equation (24a) in the main text. 




Table VII. Constants for Diffusion Coefficient* Curve Fits 



Footnotes at end of table, page 60. 




Tabic VII. Continued 


■ 

I 


Pair 

lmbcr^ 

Interaction 
pair {i - j) 

A n tJ 

% 

U IJ 

%. 

U IJ 

D n 

U IJ 

21 

NO+- NO+ 

0 

0 

3.5000 

-30.3210 

22 

e _ — N 2 

-.1147 

2.8945 

-23.0085 

65.9815 

23 

e _ — 0 2 

-.0241 

-.0029 

.3464 

.0856 

.1136 

.6655 

-1.3848 

-.8205 

24 

e” - N 

0 

0 

1.5000 

-2.9987 

25 

e~— 0 

.0581 

-1.5975 

15.4508 

-40.7370 

26 

e“ — NO 

.2202 

.2871 

-5.2261 

-8.3759 

42.0630 

82.8802 

-106.0937 

-267.0227 


27 

e — NO + 

28 

e - — e“ 

29 

N + - N 2 

30 

N+- 0 2 

31 

N+- N 

32 

N+- 0 

33 

N+- NO 

34 

N + - NO + 

35 

N + — c~ 

36 

N+- N + 

37 

0 + - n 2 

38 

0 
+ 

1 

o 

to 

39 

0 + - N 

40 

0 

+ 

1 

o 


3.5000 

3.5000 

1.9000 

1.9000 

1.5572 

1.9000 

1.8999 
3.5000 
3.5000 
3.5000 

1.9000 
1.9000 
1.9000 
1.5572 


-25.2128 

-24.8662 

-13.1144 

-13.1357 

-11.1616 

-13.0028 

-13.1254 

-30.0951 

-25.2128 

-29.9401 

-13.1578 

-13.1810 

-13.0028 

-11.1729 


Temperature 
range, K$ 


1000 to 9000 
9000 to 30 0000 


1000 to 8000 
8000 to 30000 


Footnotes at end of table, page 60. 




Table VII. Continued 



Footnotes at end of table, page 60. 



Table VII. Concluded 


Pair 

number^ 

Interaction 
pair (i- j) 

a d 

b d 

C D 

u i3 

% 

U '3 

Temperature 

range, 

60 

Oj- NO 

c 


c 


1.9000 

-13.3343 


61 

Oj- NO + 





3.5000 

-30.3036 


62 

0+— e- 







-25.2128 


63 

0+— N+ 







-30.0839 


64 

0 

1 

o 

+ 







-30.1273 


65 

0+- N+ 







-30.2867 


66 

0 
*>+ 

1 

o 

> 


' 




-30.2867 



^Diffusion coefficients are obtained in cm 2 -atm/sec. Diffusion coefficients obtained from these curve 
fits are for the limiting electron pressure p e m (eq. (23g)). For different electron pressures, the cross sections 
should be corrected by the formula given in the main text when the interacting pair of species are both 
ions or electrons or a combination of the two. Diffusion coefficients for Nj and are taken to be the 
same. 

^ Cross sections 1 to 15 are used in a 5-species air model and 1 to 28 in a 7-species model. 

tThe temperature range for all curve fits is 1000 < T < 30 000 K, except where noted. 

^Coefficient for diffusion of internal excitation energy (see p. 20). 
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Table VIII. Curve-Fit Constants for Collision Cross-Section ft-' * 



Footnotes at end of table, page 64. 










Table VIII. Concluded 



* Cross sections are obtained in A 2 , 1 A 2 = 10 _16 cm 2 . Collision cross sections obtained from these 
curve fits are for the limiting electron pressure p em (eq. (23g)). For different electron pressures, the cross 
sections should be corrected by the formula given in the main text when the interacting pair of species 
are both ions or electrons or a combination of the two. Cross sections for and are taken to be the 
same. 

t Cross sections 1 to 15 are used in a 5-species air model and 1 to 28 in a 7-species model. 

$ The temperature range for all curve fits is 1000 <T < 30000 K, except where noted. 
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— ( 2 , 2 ) 

Table IX. Curve-Fit Constants for Collision Cross-Section 


Pair 

number^ 

Interaction 
pair (i - j ) 

^ 2 ’ 2) 

V 

5 o< 2 ’ 2) 

c n (2 ' 2) 

U 

D o (2 ’ 2 > 

U ij 

Temperature 
range, ^ 

1 

n 2 - n 2 

0 


-0.0203 

0.0683 

4.0900 


2 

02 “ N2 



-.0558 

.7590 

.8955 


3 

02 ~ 02 



-.0485 

.6475 

1.2607 


4 

N - N 2 



-.0190 

.0239 

4.1782 


5 

N - 0 2 



-.0203 

.0730 

3.8818 


6 

N - N 



-.0118 

-.0960 

4.3252 


7 

0 -N 2 



-.0169 

-.0143 

4.4195 


8 

0 - 0 2 



-.0247 

.1783 

3.2517 


9 

O N 



.0065 

-.4467 

6.0426 


10 

0-0 



-.0207 

.0780 

3.5658 


11 

NO - N 2 



-.0385 

.4226 

2.4507 


12 

0 

1 

0 

2; 



-.0522 

.7045 

1.0738 


13 

NO - N 



-.0196 

.0478 

4.0321 


14 

NO - 0 



-.0203 

.0730 

3.8818 


15 

NO NO 



-.0453 

.5624 

1.7669 


16 

NO+ - N 2 



( 


-.4000 

6.7760 


17 

NO* - 0 2 










18 

NO* - N 










19 

NO + - 0 










20 

NO* - NO 








' 


21 

NO* - NO* 


- 



-2.0000 

24.3602 


22 

- N2 

.1147 

-2.8945 

24.5080 

-67.3691 



Footnotes at end of table, page 68. 



iable IX. Continued 


Pair Interaction 

number t pair (I — j) 


e O2 


c N 


e - 0 


e NO 


c NO' 


.0164 

-.2027 

-.2202 

-.2871 


-.3467 

-.0742 


-.2431 

5.6428 

5.2265 

8.3757 


C n (2 : 2) 

u 


1.3887 

.7235 


1.1231 

-51.5646 

-40.5659 

-81.3787 

- 2.0000 


V 


-.0110 

-.2116 

1.6094 

-1.5561 

155.6091 

104.7126 

265.6292 

24.3061 


Temperature 
range, K* 




1000 to 9000 
9000 to 30 000 


1000 to 9000 
9000 to 30 000 

1000 to 8000 
8000 to 30 000 


N’ - N 2 ' 
N* 0 2 


N* N 


N* O 


N- NO 


N J -NO 


N* - e 


- 2.0000 


-.4000 


-.4000 


-.4146 


-.4000 


- 2.0000 


- 2.0000 


24.3061 


6.7760 


6.7760 


6.9078 


6.7760 


6.7760 


24.3602 


24.3061 


N 4 - N + 


CP N 2 

o + 0 2 


O 4 -N 


O 4 o 


CP NO 


- 2.0000 


-.4000 


-.4000 


-.4000 


24.3602 


6.7760 


6.7760 


6.7760 


6.7787 


6.7760 


Footnotes at end of table, page 68. 






Table IX. Concluded 


Pair 

number^ 

Interaction 
pair (i - j) 

A- 

n 

(2.2) 

B n {2 ' 2) 

C f&2) 

n 

n 

(2,2) 

y 

Temperature 

range, 

60 

0 2 + - NO 

( 




-.4000 

6.7760 


61 

0 2 + - NO + 





-2.0000 

24.3602 


62 

0 2 + - e" 







24.3061 


63 

0 2 + - N* 







24.3 

602 


64 

o 2 + 0 + 










65 

o 2 + n 2 + 





1 





66 

0 2 + 0 2 + 

' 


> 


j ' 


> 















* Cross sections are obtained in A 2 ; lA 2 = 10 _16 cm 2 . Collision cross sections obtained from these 
curve fits are for the limiting electron pressure p em (eq. (23g)). For different electron pressures, the cross 
sections should be corrected by the formula given in the main text when the interacting pair of species 
are both ions or electrons or a combination of the two. Cross sections for Nj and are taken to be the 
same. " 

f Cross sections 1 to 15 are used in a 5-species air model and 1 to 28 in a 7-species model, 
i The temperature range for all curve fits is 1000 <T< 30000 K, except where noted. 
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Table X. Curve-Fit Constants for Collision Cross Section Ratio B 


(1000 K < T < 30000 K) 



Footnotes at end of table, page 72. 










Bb* 

y 

C 



D 

.19 




.44 






meters and are valid as given for 


1 to 28 in a 7-species model. 




Chemical species In high-temperature air 

Region 

Air 

chemical 

model 

Species present 

CD 

2 species 

o 2 ,n 2 


5 species 

0 2 , N 2 , 0, N, NO 

mp 

7 species 

0 2 ,N 2 ,0,N,N0, NO+.e' 

© 

11 species 

0 2 , N 2 , 0, N, NO, 0 2 , N 2 , 0 + , N + , NO + , e' 


Regions with chemical and thermal nonequilibrium 

Region 

Aerothermal phenomenon 


Chemical and thermal equilibrium 


Chemical nonequilibrium 
with thermal equilibrium 

© 

Chemical and thermal 
nonequilibrium 


i- Vibrational energy excited 


120 


Altitude, 
km 60 



Ionization 
of molecules 
0 + 0-*-0 2 + e‘ 
N + N -► N| + e" 

+ O -► NO + + e" 


Velocity, km/sec 

Figure 1. Flight stagnation region air chemistry for a 30.5-cm radius sphere (adapted from ref. 5). 
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0 


■ 10 ££ 


ln(Keq,r) 


Total number 
density 


10 1D 



10 19 


10 20 30 40 50 

Temperature, K 

(a) For different values of total number density. 

Present curve fit 

Park's curve fit (ref. 23) 

— Computed from table II 

o Exact 


60 x 10 3 


ln(Keq,r) -10 



Temperature, K 

(b) For total number density of 10 18 particles/cm 3 . 
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Figure 2. Variation of equilibrium constant with temperature for reaction N + O ~Z NO + + e 
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(a) Diatomic oxygen. 


Curve fit to Browne data 



(b) Diatomic nitrogen. 

Figure 3. Curve fit to specific-heat values obtained by Browne (ref. 35) and comparison with Hansen s values 
(ref. 37). 
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Curve fit to Browne data 



Temperafure, K 

(c) Nitric oxide. 

Figure 3. Concluded, 
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Temperature, K 

(b) Diatomic nitrogen. 

Figure 6. Curve fit to viscosity values obtained by employing collision cross sections based on data of 
reference 40. 
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Temperature, K 


(a) Diatomic oxygen. 



Temperature, K 

(b) Diatomic nitrogen. 

Figure 7. Curve fit to frozen thermal conductivity values obtained by employing collision cross sections based 
on data of reference 40. 
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Figure 8. Viscosity of equilibrium air at 1 atm. 
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Figure 9. Frozen thermal conductivity of equilibrium air at 1 atm. 
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(a) Neutral-neutral molecular interaction: O 2 «- * C> 2 - 


Computed: 



Temperature, K 

(b) Ion-neutral molecular interaction: NO + «-» 02- 

Figure 10. Curve fit to computed values of collision integrals, collision-integral ratio, and binary diffusion 
coefficient obtained by using data of reference 40. 
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Computed: 
o rf* 1 ' 1 * 



(c) Ion-ion molecular interaction NO + > NO + (p e — Pem)- 


Computed: 

o 



(d) Electron-ionized-molecule interaction: NO + e (p e = p em ). 


Figure 10. Continued. 
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